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Abstract 

We evaluate the in-medium spectral functions for pions, nucleon and isobar reso- 
nances in a self consistent and covariant manner. The calculations are based on a re- 
cently developed formulation which leads to predictions in terms of the pion-nucleon 
scattering phase shifts and a set of Migdal parameters describing important short 
range correlation effects. We do not observe significant softening of pion modes if we 
insist on reasonable isobar resonance properties but predict a considerable broad- 
ening of the iV(1440) and 7V(1520) resonances in nuclear matter. Contrasted results 
are obtained for the s-wave iV(1535) and iV(1650) resonances which are affected 
by a nuclear environment very little. The properties of slowly moving isobar's in 
nuclear matter are found to depend very sensitively on a soft form factor in the 
ttNN vertex, which is not controlled by the ttN scattering data. 



1 Introduction 



The study of pion propagation in dense nuclear matter is of central impor- 
tance when addressing the in-medium modifications of nucleon and delta res- 
onances [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. There are strong hints from the em- 
pirical photon nucleus absorption cross section as well as from pion-nucleus 
scattering data that the low-lying nucleon and delta resonances do change 
their properties in nuclear matter substantially already at nuclear saturation 
density [15,16,17,18,19,20]. Recent data on electroproduction of isobars off he- 
lium three [21] are interpreted in terms of a repulsive in- medium mass change 
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and an increased decay width of the isobar- resonance [21]. Since most nu- 
cleon and isobar resonances have a substantial decay fraction into one nucleon 
and one pion, the in-medium resonance structure reflects to a large extent 
the medium modified propagation properties of pions. The pion self energy 
in nuclear matter is quantitatively constrained by pionic atom data [6,22,23]. 
Most exciting are the recently established states where a negatively charged 
pion is bound by a heavy nucleus in an s-wave or p-wave state [24]. It is still 
an open problem to find a quantitative and microscopic derivation in partic- 
ular for the large absorptive part in the nuclear optical potential needed in 
the phenomenological description of pionic atom data [6,22,23]. The problem 
requires a non-perturbative many-body approach, since the low-density ex- 
pansion ceases to converge rapidly at the relevant nuclear densities even when 
chiral correction terms are considered [25,26,27]. 

In this work we generalize the covariant framework of [28], which was re- 
cently proposed for the self consistent propagation of antikaons and hyperon- 
resonances in nuclear matter, to the problem of pion, nucleon- and isobar- 
resonance propagation incorporating short range correlation effects. The merit 
of our scheme is that it is formulated entirely in terms of the two-body scat- 
tering phase-shifts properly extrapolated to subthreshold energies. Given the 
empirical irN phase shifts and values for the Migdal parameters [6,11,12,13] 
describing important short range correlations the scheme is parameter free. 
We expect self consistency to lead to a broadening of the pion spectral func- 
tion which may help to establish a microscopic understanding of pionic atom 
data and also offer an improved understanding of the expected broadening of 
nucleon and isobar resonances in nuclear matter. Our present approach consti- 
tutes a significant progress as compared to previous self consistent calculations 
[8,9,10] which were based on p-wave pion-nucleon dynamics only. Moreover 
this work is the first attempt to consider the effects of the in-medium mixing 
of partial wave amplitudes. In particular the feedback effect of the in-medium 
modified s-wave iV(1520) and 7V(1650), the p-wave 7V(1440) and the d-wave 
7V(1535) resonances to the propagation properties of pions will be addressed 
in this work. Similarly the isobar resonances s-wave A(1620), p-wave A(1232) 
and A(1600) and the d-wave A(1700) resonances are considered. 



2 Self consistent nuclear pion dynamics 

We briefly recall the self consistent and covariant many-body framework in- 
troduced in [28] appropriately adjusted to pion propagation in nuclear matter. 
First we recall the vacuum on-shell pion-nucleon scattering amplitude 

(ir>(q) N(p) | T |vr*(g) N(p)) = (2vr) 4 5\q +p-q-p) 



2 



x u(p) T£ N ^ N (q, p; q, p) u(p) , 



(1) 



where <5 4 (..) guarantees energy-momentum conservation and u{p) is the nu- 
cleon isospin-doublet spinor. The vacuum scattering amplitude is decomposed 
into its isospin channels. Applying standard notation for the Paul matrices Tj 
we write 

^=i/ 2) = r v > ^=3/2)=^!-^^'. ( 2 ) 

where q,p, q,p are the initial and final pion and nucleon 4-momenta and 

w = p + q = p + q, k — \{p — q) , k = \(p - q) . (3) 

In quantum field theory the scattering amplitudes follow as the solution of the 
Bethe-Salpeter matrix equation 



r dH 

T(k, k; w) = K(k, k;w) + J —-- K(k, /; w) G(l; w) T(l, k; w) , 

G(l;w) = -iS(\w + l)D(±w-l), (4) 

in terms of the Bethe-Salpeter kernel K(k, k; w), the free space nucleon prop- 
agator S(p) = l/(p' — m N + ie) and pion propagator D(q) = l/(g 2 — m 2 +ie). 
Following [29] we neglect self energy corrections in the nucleon and pion prop- 
agators. In a chiral scheme such effects are of subleading order. The Bethe- 
Salpeter equation (4) properly implements Lorentz invariance and unitarity 
for the two-body scattering process. The generalization of (4) to a coupled- 
channel system is straightforward. 

The pion-nucleon scattering process is readily generalized from the vacuum to 
the nuclear matter case. In compact notation we write 



T = JC + JC ■ G ■ T , T = T(k,k;w,u) , Q = Q(l;w,u) , (5) 

where the in-medium scattering amplitude T(k,k;w,u) and the two-particle 
propagator Q(l;w,u) depend now on the 4- velocity characterizing the nu- 
clear matter frame. For nuclear matter moving with a velocity v one has 



v c 



u. 



yi-V 2 /c 2 ' yjl - V 2 /C 2 
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We emphasize that (5) is properly defined from a Feynman diagrammatic 
point of view even in the case where the in-medium scattering process is no 
longer well defined due to a broad pion spectral function. We consider the 
effect of an in-medium modified two-particle propagator Q 



AS(p,u) = 2irie[p-u]5(p 2 - m%) {j) + m N ) Q(k 2 F + m 2 N - (u ■ p) 2 ) , 

S(p,u)=S(p) + AS(p,u), V(q,u) = — — r, 

Q(l;w,u) = —iS(\w + l,u)V(\w — l,u) , (7) 
where the Fermi momentum hp parameterizes the nuclear density p with 



In the rest frame of the bulk with = (1,0) one recovers with (8) the 
standard result p = 2k F /(3n 2 ). In a more complete approach the effect of 
nucleonic correlation and binding on AS should be considered. This is be- 
yond the scope of this work. The pion self energy Il(g, u) is evaluated self 
consistently in terms of the in-medium scattering amplitudes 



/dp _ 
-j—^ i AS (p, u)T(\(p - q),\{p - q);p + q,u) 

+ Aiife, «) , r = \ f {1 = xl2) + \ t {I =w , ( 9 ) 

where the in-medium amplitudes, 7'^ I \k,k;w,u), are defined with respect to 
the free-space interaction kernel, K, 



f = K + K-g-f = T + T-Ag-f, AQ = Q-G. (10) 

Additional contributions induced by the in-medium modification of the inter- 
action kernel /C = K + AK, are cast into the term, AII(g, u), of (9). It is crucial 
to consider the latter term since it will introduce in particular the important 
short range correlation effects described by the Migdal parameters. This will 
be discussed in more detail below. 

With (10) the self consistent set of equations (5,7,9) is rewritten in a way that 
one may start with a set of tabulated free-space scattering amplitudes 
in terms of which self consistency is achieved. Coupled channel effects, which 
are known to be important for ttN scattering at higher energies, are included 
by assigning T, Q and K the appropriate matrix structures. The effect of 



4 



all inelastic channels can be accounted for by a renormalization leading to a 
complex scattering kernel K of the 7iN channel. Therefore we may assume 
that the amplitudes in (10) already include the dynamics of all inelastic 
coupled channels. In this work we study the consequence of the in-medium 
modified nN channel exclusively. Thus it is not necessary to make the coupled 
channel structure of the nN amplitudes explicit [28]. Ultimately it would be 
desirable to also evaluate the in-medium modification of the inelastic processes 
described for instance by the nA, pN and ojN channels. 

The scattering amplitudes can be systematically decomposed into covariant 
projectors Y^\q,q\w) with good angular momentum and parity J p = (n + 
|) T . We write 



T« (M;^) = E Y( + \q,q;w) M?\yfi , n) 

n=0 

oo 

+ J2 Y n~ ) (Q,Q^)M { I -\V~s,n), (11) 

n=0 

where w 2 = s and k — | {p — q) and k = | (p — q). The representation (11) 
implies a particular off-shell behavior of the scattering amplitude, which was 
obtained in the course of constructing a systematic on-shell reduction of the 
covariant Bethe-Salpeter equation that does not depend on the choice of meson 
and baryon interpolation fields [29]. In this work we focus on the leading J — \ 

and J = | channels with Yq^ (s-wave), Yq ^ (p-wave) and y} + ^ (p-wave), Y"/ ■* 
(d-wave). For details on the construction and the properties of these projectors 
we refer to [29,28]. 

For energies above threshold \/s > m N + m n the invariant amplitudes are 
uniquely determined by the scattering phase shifts Sj\y/s) and inelasticity 
parameters T)j\y/s), 



.8nsvl l±i (^s){e - l) 

— i ± 

pl J N s + m 2 N -m 2 ±2 m N y/s 

1 ^ _ (m N + m w ) 2 ^ ^ _ (m N -m^) 2 ^ ^ 

4 s s 

We use m N = 939 MeV and = 139 MeV throughout this work. At sub- 
threshold energies y/s < — the on-shell amplitudes are again deter- 
mined by the scattering phase shifts once crossing symmetry is invoked. In 
the region m N — m n < yfs < m N + m n an analytic continuation is required as 
for instance implied by the dispersion-integral representation, 



M< ± )(V5,J-|) = 

„2 

PwN _ 

s 



5 



oo 



dw SM) ±J (w,n) 

71 W — y/s 



KM{ ±) (^,n) = ^ ±) (^,n) + ? J 



(13) 



where the functions Dj (y/s,n) represent all left-hand cut contributions [30]. 
Since the cut structure in D^\y/s^n) is largely dominated by the s- and u- 
channel nucleon-exchange contributions, it is straightforward to analytically 
continue the on-shell amplitudes into the region m N — m n < y/s < m N + m n 
as to arrive at approximate partial-wave amplitudes specified for all energies. 
As a consequence the amplitudes would receive large and strongly energy- 
dependent contributions in the interval mjv — m n < y/s < + m^. 

We point out, however, that the formulation as presented above would lead to 
unphysical and misleading results for the pion self energy unless the depen- 
dence of the partial-wave amplitudes on q 2 =/ m 2 at subthreshold energies is 
taken into account in some way. This is important, since the u-channel contri- 
bution to the partial wave amplitudes shows an extremely strong dependence 
on q 2 . Typically for given y/s with m N — < y/s < m N + m w and q 2 < 
the contribution is small and negligible but at the on-shell point q 2 = m 2 it 
is large. Since the subthreshold amplitudes are probed by the pion self en- 
ergy only for q 2 < m 2 , it is clear that the g 2 -dependence must be considered. 
Rather than extending our scheme for the most general off-shell structure in 
the scattering amplitude we propose a simple modification of the scheme that 
successfully circumvents the artifacts of a pure on-shell scheme. It would be 
unclear in any case what to use for the off-shell dependence, since such a 
dependence is highly scheme-dependent. 

The idea the proposed scheme is based on, exploits crossing symmetry as a 
tool that determines the troublesome off-shell part of the scattering amplitude. 
We seek a representation of the free-space scattering amplitude of the form 



that manifests the desired crossing symmetry of the amplitudes explicitly but 
masks the constraints set by unitarity. The decomposition (14) is not unique 
per se but can be tailored to our application. In particular we insist that T s 
contains only s-channel unitarity cuts but shows no cuts in the u-channel. 
Similarly we insist that T u is real for y/s > uin + rn n . In general there will 
be contributions present in T which lead to both an s-channel and u-channel 
cut. Such contributions must be put into a remainder term not displayed in 
(14), if we wish to keep perfect analytic properties of the amplitudes T s and 
T u . We point out however, that such contributions are zero identically, if one 
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considered elastic scattering processes only. Therefore we expect the degree of 
analyticity violation, when demanding T s to represent the complete strength 
of the s-channel unitarity cut, to be reasonably small if energies are not ex- 
cessively large. Within a given application window there is still the freedom 
how to distribute a polynomial background contribution. Here it is advanta- 
geous to insist on T s = T beyond a certain scale \/s > A. Consistency then 
requires that for y/s < (m 2 N — m^)/A one should obtain T s = or equivalently 
T = T u . The scattering amplitude T s is decomposed into partial wave ampli- 
tudes Mj ± \^, n). In practise we use the representation (13) with the function 
^t\\^-> n ) replaced by a second order polynomial in y/s. The coefficients of 
the latter are determined to guarantee first, that the real parts of Mj^-^/s, n) 
vanish identically at \/s = (m% — m 2 )/ 'A and second, is such that a smooth 
matching to the real part of the partial wave amplitudes, Mj ± \\/s,n), at 
\/s = A as given in (12) is obtained. Of course in the Pu channel the nucleon 
pole contribution is incorporated in addition. The pion-nucleon coupling con- 
stant, g w NN = 13.34, of [31] is used. With A = 1600 MeV a remarkably smooth 
matching is achieved for all considered channels as described above. The de- 
viation of Mj^-^/s, n) from M\ ± \^,n), both shown in Fig. 2, is a measure 
for the importance of u-channel cuts in a given partial wave. 

The idea is to use the amplitudes T s in (14) properly decomposed into partial 
waves via (11). The effect of T u is incorporated into the pion self energy by 
adding to (9) the term induced by T s but with q — > —q replaced. As we 
apply the decomposition (14) for a many-body evaluation of the pion self 
energy, the parameter A plays the role of a scale that determines at what 
energy the manifestly crossing-symmetric scheme recovers two-body unitarity 
exactly. For energies close to the pion-nucleon threshold a unitarization is not 
really required for a quantitative description of the scattering amplitude. The 
amplitudes are largely dominated by the s- and u-channel baryon exchange 
contributions. In contrast, at higher energies the unitarity constraint becomes 
more and more important rendering any constraint from crossing symmetry 
rather implicit. 

The solution of the in-medium Bethe-Salpeter equation is considerably com- 
plicated by the presence of further tensor structures and Q\ij\ not required 
in the vacuum. For a complete collection of the in-medium projectors together 
with their induced in-medium loop functions, AJ^\w,u), representing the 
object AQ in (10), see [28]. In this work the loop functions are evaluated by 
their dispersion-integral representation in terms of their imaginary parts. A 
subtraction constant is determined by insisting that the scattering amplitudes 
do not renormalize the nucleon mass 



A ^fW) ,,_ 2 =0 - ry\q,q-w,u) = TP(q,q,w), (15) 
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Fig. 1. Partial wave irN scattering amplitudes at subthreshold energies as recon- 
structed from the SM01 phase shifts [32]. The solid and dashed lines give the 
real and imaginary parts of the M^ iz \y/s, J - 1/2) amplitudes of (12). The addi- 
tional solid lines extending down to subthreshold energies represent the amplitudes 
Mj^\y/s, J — 1/2) that represent Ts J \q, q; w) of (14) and are free of u-channel cuts. 

for w ■ u > 0. Such a condition is well justified in the present scheme since 
neither correlations nor binding effects are incorporated. That is outside the 
scope of this work. Furthermore we drop a small contribution of intermediate 
hole states, not consistently treated here in any case, but part of the relativistic 
scattering equation (5). This leads to an in-medium scattering amplitude, 
whose imaginary part, 



^sT s (q, q;w,u) = if w ■ u < \j m 2 N + k 2 F , (16) 

vanishes below the chemical potential as expected from the Pauli principle. 
To be precise the property (16) requires also the consideration of the Pauli 
blocking effect in the s-channel nucleon exchange contribution, i.e. IC ^ K. 

We turn to the effects of short range correlation [6] . Here we follow the recent 
work [14] and apply covariant expressions parameterized in terms of the Migdal 
parameters g'u,g'i 2 an d #22 • The delta-hole term is folded by an isobar spec- 
tral function obtained from the appropriate in-medium scattering amplitude 
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assuming the averaged value w = 200 MeV. We then approximate the isobar 
spectral function to be a function of Wq—w 2 . Thus, the value for g' 22 used in this 
work is introduced with respect to the in-medium pion-isobar coupling con- 
stant. The latter will show a sizeable in-medium reduction in our scheme. The 
expression for the pion self energy as given in [14] properly subtracted by its 
leading order contribution, defined in the limit g'^ — > 0, is identified with All 
in (9). The real part of this contribution is evaluated in terms of a dispersion- 
integral representation with a subtraction constant fixed at oj 2 — q 2 = m 2 
matched to the real part of the original All contribution. Since the nucleon- 
hole contribution as derived in terms of the in-medium scattering amplitude 
includes an appreciable momentum dependent renormalization of the pion- 
nucleon coupling constant, a corresponding renormalization is applied in All. 
This effect reduces the strength of the nucleon-hole contribution. The values 
we quote here for Migdal's g' u and g' 12 parameters are defined with respect 
to the free-space pion-nucleon but in-medium pion-isobar coupling constants 
[14]. 

We reconstruct the real part of the pion self energy, $tll(q, u), in terms of 
its imaginary part, SII(g, u). From now on we work in the rest frame of nu- 
clear matter with = (1,0) for convenience and write u) = U(u,q). A 
subtracted dispersion-integral representation is imposed 



t-t/ -.x ,^ . } d^ 2 3n(w,g) 

U{uj,q) = c{q)+ — -, (17) 

J ix uj z — ur — i e 



cut off by A = 1.2 GeV. The reflection property II (u;, q) = II(— u, q) holds for 
isospin symmetric nuclear matter as a direct consequence of our manifestly 
crossing-symmetric scheme. The subtraction constant c(q) is fixed such that 
the real part of the self energy reproduces the corresponding value for the 
contribution of the amplitude defined by (9) at the point 

oj = Max^AS + dgl + M 2 " \A< + k h \f^ 2 + W\ 2 ~ ™ N ] . 

This condition guarantees that the real part of the in-medium amplitude is 
probed only for w 2 , — w 2 > A 2 , where the effective amplitude T^(q,q;w) 
represents the exact scattering amplitude (see (def-decom)). We assure that 
the sum rule 



/ — 3 D(u, q) = — 1 , 

J IT 



(18) 



holds accurately within 1 % in our scheme. 
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The self consistent set of equations (9,5) is solved numerically by iteration 
where we start with the evaluation of the in-medium self energy. Convergence 
is typically found after 5 to 6 iterations. The energies and the three-momenta 
are restricted by \lu\ < 1.2 GeV, \q \ < 0.8 GeV and w < 2 GeV, \w\ < 1.2 
GeV. The free-space partial wave amplitudes, Mj ± \^, n), are put to zero for 
y/s > 2 GeV and yfe < 0.54 GeV. 



3 Results 

We give a presentation and discussion of our results for the in-medium modi- 
fication of the pion and the J = |, | nucleon and isobar resonance properties. 
The resonance propagator can be identified with the appropriate pion-nucleon 
scattering amplitude of a given partial wave. In the self consistent scheme of 
section 2 the in-medium scattering process is intimately related to the pion 
spectral function. According to (9) once the self consistent in-medium scat- 
tering process is established the pion self energy follows by averaging the 
in-medium scattering amplitudes (9) over the Fermi distribution. Therefore 
the pertinent structures in the in-medium amplitudes already tell the charac- 
teristic features expected in the pion spectral function. 

Since the present day literature does not offer a unique set of Migdal parame- 
ters g'u, g 22 and g' 12 we studied first the sensitivity of our approach to different 
choices thereof. We use two sets of parameters in this work. Set I) is assuming 
9n = 9i2 = 922 = 0-8- Set II) is given by the values g' n = 0.585 , g' 12 = 0.15 to- 
gether with g' 22 = 0.60. We remind the reader that we use here a convention in 
which g' n is defined with respect to the free-space pion-nucleon coupling con- 
stant whereas g' 22 is defined with respect to the in-medium pion-nucleon-isobar 
coupling constant. 

With the nuclear density set to p = 0.17 fm~ 3 Fig. 2 shows the resulting pion 
self energy together with the real and imaginary parts of the in-medium ttN 
scattering amplitude in the P 33 channel. As illustrated by the figure the ef- 
fect of different choices for the Migdal parameters affects the pion self energy 
mostly at not too large energies and intermediate momenta. As illustrated 
by the right hand panels the effects of such differences are quite important 
for the properties of the A(1232)-isobar in nuclear matter. Only set I) leads 
to acceptable properties of the isobar though we may still overestimate the 
broadening of its width in nuclear matter for high velocities. Note that the 
depletion of the isobar peak in Fig. 2 is a combined effect of reduced pion- 
isobar coupling constant and increased decay width. Thus the g 22 parameter 
as defined more conventionally with respect to the free-space pion-nucleon- 
isobar coupling constant is in fact smaller than 0.8 about 0.6 only. For isobar 
propagation with non-zero momenta set II) leads to a substantial broaden- 
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Fig. 2. In-medium pion self energy (left hand panels) and ttN scattering amplitude 
in the P33 channel (right hand panels) at nuclear saturation density. The solid lines 
show the results obtained with set I) of Migdal parameters, the dashed lines those 
for set II) as specified in the text. The upper two (lower two) right hand panels 
describe isobar propagation in nuclear matter with total isobar three momentum 
\w\ = 0.0 GeV (\w\ = 0.4 GeV), compared to free space results. 

ing unlikely to be consistent with photo-absorbtion data of the nucleus [16]. 
We confirm the results of [3] that the splitting of transverse (Q-space) and 
longitudinal (P-space) modes is small. In our formulation such effects reflect 
the in-medium mixing of partial wave amplitudes. In none of the channels we 
observed significant effects thereof. The attractive mass shift for the isobar of 
about 80 MeV is in qualitative agreement with previous microscopic calcula- 
tions [3]. Referring to the argument put forward in [1,3] the apparent repulsive 
mass shift as extracted from photo absorption data [4] is a combined effect of 
an attractive isobar self energy and short range correlation effects. 

The results differ from conclusions of previous works [8,9,10] that also incor- 
porated self consistency but did not observe such a significant broadening of 
the isobar. This difference is a direct consequence of the soft form factor in the 
pion-isobar vertex used in [8,9,10]. In our scheme there is not much place for 
such phenomenology, since the applied phase shifts entail already form factor 
effects that one may want to use modeling the pion-nucleon interaction. We 
emphasize our goal to develop a microscopic understanding of pion propaga- 
tion in nuclear matter. First we continue to present and discuss results solely 
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based on the ttN phase shifts and a set of Migdal parameters only. Never- 
theless, we will return to this issue and also present results that follow upon 
incorporating a phenomenological form factor into our scheme. 



3. 1 Propagation without form factors 

In the form factor-free scheme the isobar width can be protected against ex- 
cessive increase to some extent by an appropriate choice of Migdal parameters. 
Large g' 12 and g' 22 parameters together with a reasonably large g' u ~ 0.6 — 0.8 
lead to a suppression of fast soft modes that are responsible for the broad- 
ening of the isobar, i.e. the low-energy tail of the pion spectral function at 
momenta \q\ ~ 400 — 500 MeV is suppressed. We are aware that this con- 
clusion is possibly inconsistent with the conclusion of the recent work [33]. 
This issue certainly deserves more detailed investigations. Vertex correction 
diagrams that mask the pion spectral function as probed in the isobar self 
energy but not considered in the present scheme yet may lead to a similar 
suppression of fast soft modes. Thus we would not exclude a small value of 
g' 12 once such effects are incorporated. 

A further interesting question is the amount of softening found in the pion 
self energy. Typically the traditional nucleon and isobar-hole model leads to 
a minimum of the function 



at some intermediate momentum. Such a minimum helps to explain for in- 
stance unnatural parity states of finite nuclei [6]. To be specific the covariant 
nucleon and delta-hole model proposed in [14] shows a minimum at \q\ ~ 365 
MeV and nuclear saturation density using the parameter set II). In the present 
self consistent scheme we do not find a minimum of S(q) for either of the two 
choices of Migdal parameters studied here. The size of the function at a given 
momentum, say \q\ = 300 MeV, strongly depends on the choice of Migdal 
parameters, with S ~ 3.4 m\ and S ~ 2.4 m\ for set I) and II) respectively. 
By lowering g' u and g' l2 down to unrealistic values it is possible to produce a 
minimum of S(q) in our present scheme, however at the prize that the isobar 
resonance is dissolved almost completely. 

At small pion momenta \q\ < 40 MeV the pion self energy may be approxi- 
mated by 



S(q) =q 2 + ml + U(0,q), 



(19) 



• 2 =5R[am2 +/ 3 5*2] 



(a-l)ml + (P-l)q 



(20) 
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Fig. 3. Real (left hand panel) and imaginary (right hand panel) part of the pion self 
energy at p = 0.5 po, po and 1.5 po with po = 0.17 fm" 3 . The results were obtained 
with parameter set I). 

the values of the parameters known to some extent from pionic atom data [5,7]. 
For set I) and II) we obtain (a, (3) ~ (1.09 - i 0.05, 0.56 - i 0.25) and (a, (3) ~ 
(1.06—i 0.07, 0.74— i 0.44). We do not find any significant effect from a possible 
wave function renormalization. Comparing (20) with the results of [7] the s- 
wave absorption strength parameterized by S a comes out about a factor two 
too small. This reflects the fact that not all absorption channels are included 
in our work. The value for 9ft a is reasonably close to phenomenological values. 
The p-wave absorption strength related to 9 f3 is overestimated somewhat [5]. 
This may be linked to a possibly too large in-medium isobar width obtained 
in our present scheme. The parameter 9ft /3 of set I) is however close to the 
empirical value [5]. Given the fact that our scheme is parameter free except for 
the choices of the Migdal parameters, these results are encouraging. However, 
it should be emphasized that one expects some non-linear dependence of a and 
j3 on the density. For instance at 0.5 p we obtain (a, (3) ~ (1.09 — i 0.02, 0.46 — 
i 0.11) for set I). To further illustrate the amount of non-linear behavior in the 
pion self energy we include Fig. 3 which displays the results for the pion self 
energy for the nuclear densities 0.5 po, 1.0 po and 1.5 po- A striking non-linear 
behavior is typically seen at small pion energies. 

In Fig. 4 our results for the pion spectral function and the J = | isobar 
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Fig. 4. Pion spectral function (left hand panels) and J = | nucleon resonance prop- 
agators (right hand panels) at po = 0.17 fm -3 (solid lines) and 2 po (dashed lines) as 
functions of the pion, u, and resonance, wo, energy respectively. The vertical lines in 
the left hand panels show the energy of a pion in free space at given momentum \q\. 
The upper (lower) two right panels show the in-medium ttN scattering amplitude 
in the P33 (-D13) channel. Results are shown for a resonance at rest with \w\ = 
(thick lines) and for three-momentum \w\ = 0.4 GeV (thin lines) with longitudinal 
polarization. 

A(1232) and iV(1520) resonances are shown for set I) of Migdal parameters 
at two different nuclear densities, po = 0.17 fm~ 3 and 2 p . Of course the 
results at 2 p should be considered cautiously because nuclear binding and 
correlation effects were not yet fully included in the present scheme. The pion 
spectral function clearly exhibits the three well known modes, zero-sound, 
isobar-hole and pion branch with weighting factors strongly dependent on the 
pion momentum. We wish to make two points here. First, in contrast to the 
standard delta-hole model for the pion [1] , density we do not observe a signifi- 
cant strength of soft pion modes at nuclear saturation in our work. This is due 
to self consistency. Using somewhat smaller values gL = 0.6 does not change 
our conclusion qualitatively. The resulting spectral function is quite similar to 
the one shown in Fig. 4 only that for instance at \q\ = 400 MeV about 25 % 
of the strength sitting in the pion branch is moved up to the delta-hole branch 
and to a lesser degree into the zero-sound branch. Our finding confirms the re- 
sults of previous self consistent approaches to pions in nuclear matter [9,8,10] 
qualitatively and should have important consequences in various applications 
of the pion spectral function to hadron properties in nuclear matter. Second, 
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the significant broadening and repulsive shift of the main mode at large mo- 
menta is obviously a result of the inclusion of s- p- and d-waves in our scheme. 
Such effects are lacking in a model that incorporates nucleon- and delta-hole 
terms only. 

We turn to the resonance properties. Fig. 4 shows the narrowing of the iso- 
bar resonance with increasing density due to the growing importance of Pauli 
blocking. A similar but much less pronounced effect is observed for the iV(1520) 
resonance. The significant broadening of the latter resonance at nuclear satura- 
tion density is in qualitative agreement with previous phenomenological works 
[17,18,19,20] and may help to arrive at a microscopic understanding of photo 
absorption data in the second resonance region expected to be dominated by 
the contribution of the iV(1520) resonance [34,35,36]. The fact that at nuclear 
saturation density we do not obtain any large in-medium effect from the nN 
channel on the remaining resonances, JV(1535), iV(1650), A(1620), A(1600) 
and A(1700) except for the iV(1440), for which we find a significantly increased 
width, is interesting and deserves further studies. The Di 3 and Pu resonances 
iV(1520) and iV(1440) appear particularly sensitive to the in-medium dress- 
ing of the nN channel due to their d-wave and p-wave phase space behavior. 
Clearly the effect of possible in-medium modifications of inelastic channels like 
7rA, pN, ujN on the resonance properties asks for further detailed studies [37]. 



3.2 Propagation with form factors 



Until now we did not include any phenomenological form factors into the 
computation, motivated by the fact that most of such effects are already taken 
care of by using the physical scattering amplitudes. However, a certain form 
of pion-nucleon-nucleon and pion-nucleon-isobar form factors commonly used 
in calculations of pion, nucleon and isobar self energies based on the relevant 
3-point functions can also be introduced in the present scheme. The condition 
is that the solution of the Bethe-Salpeter equation in vacuum is not affected 
for the on-shell scattering amplitudes. Indeed, the use of scattering phase 
shifts incorporates the form factors necessary to describe the pion-nucleon 
interaction, but only for on-shell particles. Off-shell-pion effects may allow the 
inclusion of phenomenological form factors of the form, 

PnNx(q 2 ) = exp \-{q 2 - m^f/A^ J , (21) 



which suppress off-shell-pion contributions from relevant channels (with X 
being N or A) when calculating the pion self energy. It should be emphasized, 
however, that the phenomenological need for such a form factor may stem 
from medium modifications of the interaction vertices rather than a strong 
off-shell dependence of the vertices in free-space. We do not modify the loop- 
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Fig. 5. In-medium ttN scattering amplitude in the P33 channel at nuclear saturation 
density, compared to the vacuum amplitude (solid line) . The values of Migdal param- 
eters are: g' n = g' 12 = g'22 = 0-8. The dashed line corresponds to A w nn = 0.5 GeV 
and Ajttva — ► 00, the dash-dot line to A^nn = A w na = 0.5 GeV, the dash-dot-dot 
line to the calculation without form factors. The left two (right two) panels describe 
isobar propagation in nuclear matter with total isobar three momentum \w\ = 
GeV (\w\ = 0.4 GeV). 

integral calculation, since for consistency that would require considering also 
the influence of the Migdal parameters on the pion-nucleon-nucleon and pion- 
nucleon-isobar vertices in those loop integrals. This is not considered in this 
work due to serious complications. 

In Fig. 5 we show the effect of the form factors on the irN scattering am- 
plitude in the P 33 channel at nuclear saturation density (other channels are 
affected marginally). Whereas the form-factor influence is quite pronounced at 
momentum |iZ? j = it is of minor importance at \w\ = 0.4 GeV. As expected, 
the nNN form factor reduces the strength of the pion spectral function at 
low energy (less than m n ), thus suppressing the isobar decay amplitude and 
decreasing its width. The influence of the 7riVA form factor is less pronounced, 
since the delta-hole branch of the pion spectral function is not so far off-shell. 
In the present approach it somewhat broadens the isobar and also makes its 
mass shift less pronounced. 

The pion spectral function is not dramatically changed by the form factors, 
whose effect is illustrated in Fig. 6. The figure shows the effect of the 7rNN 
form factor since it produces the most pronounced effect. In general we get 
slightly more softening of the main pion mode and a suppression of the very- 
low-energy part. The pion- nucleon- isobar form factor does not bring in no- 
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Fig. 6. In-medium pion spectral function at nuclear saturation density, for different 
values of pion momentum: (solid line), 200 MeV (long dash), 400 MeV (short 
dash), 600 MeV (dot-dash line). The Migdal parameters are g' n = g[ 2 = g'22 = 0-8- 
The black line corresponds to calculation without form factors, while the gray lines 
are obtained with A w nn = 0.5 GeV, A. n NA — ► 00. For zero momentum the change 
is marginal and not shown. 



ticeable softening, although it makes the medium-energy part of the spectral 
function more filled in, without suppressing the small-energy region. 

In Fig. 7 we show the pion self energy at nuclear saturation density, for mo- 
menta \q\ = 0.2 GeV and \q\ = 0.4 GeV. The most prominent feature of the 
self energy is the strong suppression of the delta-hole contribution for mo- 
menta around 0.2 GeV, as well as the contributions of higher-mass resonances 
which leave the imaginary part non-vanishing even at energies up to 1 GeV. 

These results show that further investigations of the effects of vertex modifi- 
cations in the medium are required, which respect the features intrinsic to our 
present scheme. 



17 




Energy [GeV] Energy [GeV] 



Fig. 7. In- medium pion self energy at nuclear saturation density, for |g| = 0.2 GeV 
and \q\ = 0.4 GeV pion momentum. The full line corresponds to a calculation 
without form factors with g' u = g' 12 = g'22 = 0-8. The dash-dot line shows the 
results for A w nn = A-kNA = 0.5 GeV and unchanged g' values. For comparison we 
show the results (dashed line) of a non-relativistic computation from ref . [9] . 

4 Summary 

In this work we evaluated the pion self energy in nuclear matter in a self 
consistent and covariant manner. A novel framework based on the 71N phase 
shifts, as measured in free space, together with a set of Migdal parameters was 
developed. Important constraints of crossing symmetry and unitarity were in- 
corporated approximatively. Using reasonable values for the Migdal parame- 
ters we found that the nucleon resonances A(1535) and iV(1650) are basically 
unaffected by the nuclear environment. Contrasted results were obtained for 
the p-wave A(1440) and d-wave 7V(1520) resonances for which we predict 
considerable broadening already at nuclear saturation density. 

Our result for the isobar resonance are not satisfactory at this stage, due 
to a significant overestimate of its in-medium decay width. Improved results 
were obtained by incorporating a soft phenomenological form factor into the 
irNN vertex. Whereas the properties of slow isobars in nuclear matter are 
changed significantly, a soft form factor has rather moderate effects on the 
pion spectral function. These results show that further detailed investigations 
of the effects of vertex modifications in the medium are required to arrive at a 
fully microscopic understanding of the properties of isobars in nuclear matter. 
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